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^ ■ Abstract 

. Theoretical calculations for the VF-pair production process at LEP2 in terms of 

i Monte Carlo event generators RacoonWW and KorlaW&YFSWWS are reviewed. The dis- 

' cussion concentrates on precision predictions for the main LEP2 WW observables. 

r"| ' The theoretical precision of the above programs is estimated to be ~ 0.5% for the 

Q-i! total WW cross section aww, ^ 5MeV for the VF-boson mass M\y, and ~ 0.005 

Q^l for the triple-gauge-boson coupling A = = A^, which is sufficient for the final 

^ ■ LEP2 data analyses. 
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1 Introduction 



The process of VT-pair production in electron-positron colliders is very important for 
testing the Standard Model (SM) and searching for signals of possible "new physics" ; see 
e.g. Refs. . One of the main goals when investigating this process at present and future 
e^e~ experiments is to measure precisely the basic properties of the W boson, such as its 
mass My/ and width Vw- This process also allows for a study of the triple and quartic 
gauge-boson couplings at the tree level, where small deviations from the subtle SM gauge 
cancellations can lead to significant effects on physical observables - these can be signals of 
"new physics" . On the other hand, the process ly-pair production and decay constitutes 
a significant background to other processes, in particular to Higgs-boson production 

During the 1996-2000 LEP run, the four LEP2 experiments collected about 10,000 
VF-pair events each. This allows a test of the VT-boson physics at the precision level 
of ~ 1% @]. The main WW observables investigated at LEP2 are (1) the total WW 
cross section awwi (2) a distribution of the W invariant mass, which is used to extract 
the H^-boson mass Mw, and (3) a distribution of the W production polar angle cosOvf, 
which is the most sensitive observable to the triple-gauge-boson couplings (TGCs). The 
expected experimental accuracies of the final LEP2 data analyses for aww^ Mw and the 
C- and P-conserving TGC A = = are given in Table 1. In this table we also present 
requirements for the theoretical precision of these quantities. For the theoretical errors 
(Sth) it is commonly required that they should be smaller than the experimental ones (Sex) 
by a factor of at least 2. 

Table 1: 



Quantity 


Final LEP2 6ex 


Required 6th 


(Tww 

Mw 
A 


~ 1% 
- 30 MeV 
~ 0.01 


< 0.5% 
< 15 MeV 

< 0.005 



The final LEP2 data analyses need theoretical predictions in terms of Monte Carlo event 
generator (s) (MCEG) that meet these precision requirements. 

2 Theoretical description of VF-pair production 

The basic lowest-order process in which W pairs were produced at LEP2 is: 

e+ + e- -.W+ + W- -^fi + f2 + f3 + U (1) 

This can be described by three Feynman diagrams, called CC03: two annihilation dia- 
grams with 7 and Z in the s-channel (this is where the TGCs appear), and a conversion 
diagram with an electron neutrino in the t-channel (this contribution dominates at LEP2 
energies). And here the first theoretical problem appears. Namely, these three diagrams 
alone are not gauge-invariant. In order to obtain a gauge- invariant scattering amplitude, 
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one has to include a certain class of singly l^-resonant diagrams. Therefore, the min- 
imal gauge-invariant subset of Feynman diagrams containing W pairs is the so-called 
CCll-class - for more details see e.g. 0. This leads in practice to the necessity of consid- 
ering the full 4/-process, i.e. e~^e~ — 4/, which considerably complicates the theoretical 
description. 

The next problem that appears in the theoretical description of this process is the 
inclusion of a finite VT-boson width Tw- This inclusion is necessary to avoid singularities 
in the scattering amplitude coming from 14^-boson propagators. Up to now there is no fully 
satisfactory approach to deal with this in the general 4/-process. Some simple schemes, 
such as the "fixed-width scheme" and the "running-width scheme" are known to violate 
gauge invariance. The schemes that do not violate gauge invariance, such as the "fermion- 
loop scheme" or the "complex- mass scheme" possess other drawbacks, see e.g. In 
practice one commonly uses the "fixed-width scheme", for which gauge- violating effects 
are numerically negligible at LEP2 energies 0. 

All the above problems increase dramatically when one goes to 0{a) calculations. At 
the one-loop level, the number of Feynman diagrams increases to several thousands per 
4/-channel (final state) and the inclusion of the finite py-boson width gets much more 
difficult. 

One may ask the question: Do we need to go beyond the Born level? So let us briefly 
discuss the various radiative corrections and their effects on the main LEP2 WW mea- 
surements. 

Pure QED correction: 

QED corrections can be divided into a few classes: (1) initial-state radiation (ISR), i.e. 
photon radiation from incoming beams; (2) Coulomb correction - electromagnetic inter- 
action of slowly moving W^W~] (3) final-state radiation (FSR), i.e. photon radiation in 
W decays; and (4) non-factorizable (NF) corrections corresponding to interconnections 
of various stages of the process. They affect the WW observables in different ways. The 
most sizeable numerical effects come from the ISR. For aww they amount to between 
about —20% near the WW threshold and about —5% at -Ecm = 200 GeV (they depend 
stronlgy on Ecm)- Mw is affected at ~ lOMeV, but this may be enhanced by experi- 
mental reconstruction effects (kinematic fits, etc.). The TGC A is shifted by the ISR by 
~ 0.07. The Coulomb correction modifies mainly aww at the level of ~ 6% near the 
WW threshold. The FSR infiuences considerably the Mw measurement, at the level of 
~ 10-80 MeV, which strongly depends on experimental acceptances The NF correc- 
tions affect only Mw at the ~ 1-5 MeV level, when treated inclusively (which is a good 
approximation for LEP2 experimental acceptances). 

Electroweak (EW) corrections: 

The leading EW corrections connected with an effective scale of the H^-pair production 
and decay process can be taken into account by using the so-called scheme. This 
corresponds to parametrizing the cross section by the Fermi constant instead of the 
fine-structure constant a. Numerically, this changes the overall normalization, i.e. aww 
by ~ 15%. Then the 0{a) non-leading (NL) EW corrections amount to 1-2% for aww 
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at LEP2 energies. They also affect the TGCs, e.g. A at the level of ~ 0.01-0.02. 
QCD corrections: 

QCD corrections affect normalizations as well as events shapes of the hadronic WW 
channels. They are usually accounted for by including the so-called naive QCD correc- 
tion, i.e. the normalization factor as(M^)/7r ~ 3.8% for each final-state quark pair and 
2as{M^r)/37^ for Tw, in the parton-level calculation, while leaving the exclusive QCD 
effects to be modelled by dedicated MC packages, such as PYTHIA, etc. 

All the above effects are necessary in a MCEG for l^-pair production at LEP2. The 
most difficult are 0{a) NL EW corrections. To date, there exist no complete one-loop 
calculations, even for the simplest CCll-type channels. Even if such calculations existed, 
they would probably be very complex and slow in numerical computation. Therefore they 
would not be useful for MC event generation. These are the reasons why some efficient 
approximations for the WW process are necessary, at least for LEP2. 

In the case of ly-pair production, one is interested in the doubly ly-resonant process. 
Therefore one can make use of another expansion parameter, which is Tw/Mw ~ 2.5%, 
and apply the so-called pole expansion, i.e. an expansion about the complex pole cor- 
responding to an unstable W. This is the expansion in increasing powers of Tw/Mw, 
which corresponds to decreasing powers of the resonance. Then, one usually applies the 
so-called leading-pole approximation (LPA) in which only the highest-pole (resonant) con- 
tributions are retained, i.e. terms ~ (Tiy/M^y)^. For the VT-pair production process this 
means retaining only double-pole contributions, i.e. applying a double-pole approxima- 
tion (DPA). In this context LPA means just DPA; however, LPA itself has a more general 
meaning. The resulting matrix element is gauge- invariant, and the imaginary part of 
the pole position corresponds to the finite IV-boson width, see e.g. In general, 

the LPA is however not sufficient at the Born level, because the error introduced by this 
approximation can be here ~ Tw/Mw, i.e. ~ 2.5%. It may also not be sufficient for the 
leading corrections because of big-log enhancements, i.e. ~ (T\y/Mw){c(/7i) log(Q^/m^), 
where Q is some large momentum transfer and m is a small fermion mass. Here lower-pole 
terms might be necessary for the desired experimental precision. In practice, it is usu- 
ally simpler to take the whole e'^e~ —> 4/, including the leading corrections. Then, 
for the genuine 0{a) NL EW corrections the uncertainty introduced by the LPA is 
~ (Tw/M\Y){a/7i) < 10~^, therefore negligible for the LEP2 precision. 



The above solutions have been implemented in two MC event generators: RacoonWW 10- 



1^ and KoralW&YFSWWS [[7|,p!^-p2|. These programs are briefly described in the next two 



sections. 



3 The MC program RacoonWW 

The main features of the program RacoonWW are: 

• Matrix elements for all e~^e~ 4/ and e'^e^ — > 4/7 processes in massless-fermion 
approximation. 
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• The ISR in the leading-log (LL) approximation with soft-photon exponenti- 
ation through QED structure functions. 

• Coulomb correction for off-shell W^s. 

• The NF virtual corrections in the DPA. 

• The virtual 0{a) NL EW corrections in the DPA (for which one- loop calculations 
of on-shell WW production and decay are used). 

• Two QCD-inspired methods of treating soft and collinear photon singularities: dipole 
subtraction and phase-space slicing (to get a proper matching between virtual and 
real-photon corrections) . 

• Anomalous Triple and Quartic gauge-boson couplings. 

• Multichannel MC algorithm for integration and event generation. 



More details can be found in Refs. 10-13 



4 The MC programs KoralW and YFSWW3 

KoralW and YFSWW3 are actually two independent MCEGs, with some specific and some 
common features, which are listed below. 

KoralW specific features: 

• The fully massive matrix element for all e~^e~ —>■ 4/ processes (generated by the 
GRACE System |23|). 

• Two independent efficient multichannel MC algorithms for the full 4/ phase space. 

YFSWW3 specific features: 

• Multiphoton radiation in the WW production stage in the Yennie-Frautschi-Suura 
(YES) exponentiation scheme. 

• The 0{a) NL EW corrections in the LPA (based on the 0{a) calculations for on- 



shell WW production of Refs. P4P5[]). 
Common features: 



• The ISR in the YES exponentiation scheme up to 0{a^) LL with non-zero pt multi- 



photon radiation. 

Coulomb correction for off-shell W^s. 



• Non-factorizable QED corrections in an inclusive approximation of the "screened- 
Coulomb" ansatz p. 

• Anomalous TGGs (three parametrizations). 
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• The FSR at C(a^) LL generated by PHOTOS |^ (with non-zero pt photons). 



r decays done by TAUOLA quark fragmentation and hadronization managed by 
JETSET [|8|. 



• The semi-analytical program KorWan for the WW process including leading correc- 
tions (for test, fits, etc.). 

As can be seen from the above description, KoralW is devoted to the full e~^e~ — > 4/ 
process and it includes the leading corrections, i.e. all except 0{a) NL. YFSWW3, on 
the other hand, is dedicated to a precision description of the signal WW process, i.e. it 
includes all numerically important radiative corrections to this process treated in the LPA, 
but does not include the 4/-background contributions. However, thanks to the common 
features of the two programs, it is easy to correct the predictions of one of them, so that the 
final results include all the necessary contributions. In fact, the two programs implement 
special reweighting routines for this purpose. Moreover, we managed to combine the 
results of the two programs in a real-time execution (event-by-event) using the Unix/Linux 
named (FIFO) pipe mechanism. This resulted in the so-called Concurrent Monte Carlo 
(CMC) KoralW&YFSWWS. This CMC works effectively as a single MCEG and provides the 
description of the W^-pair production process with all contributions/corrections needed 
for the LEP2 precision - for more details see Ref. 



5 Theoretical precision of the main LEP2 WW ob- 
servables 



The precision requirements of the theoretical predictions for some of the main LEP2 
observables has been shown in Table 1. Any MCEG to be used for the final LEP2 
data analysis should satisfy these requirements. Estimating the precision of theoretical 
calculations is a difficult task, because it requires assessing the size of missing higher-order 
contributions. In practice, this is done by comparing different calculations (programs), 
switching on and off various contributions/corrections and looking at the effects for a given 
observable. In some cases this can be supplemented with a simplified scale-parameter 
analysis. All such investigations can give a hint on the size of the missing effects in the 
theoretical description. 

The programs RacoonWW and YFSWW3 have been extensively cross-checked and com- 
pared for many observables, see e.g. [^]. It has been found that for aww the two pro- 
grams agree within 0.3% over the LEP2 energy range and they describe the LEP2 data 
very well |^ . From the comparisons between these two programs, the comparisons with 



other calculations and investigations of various effects, the theoretical precision for aww 
at LEP2 has been estimated at ~ 0.5% Similar studies, although more involved, 
were also performed for the W mass and the TGC A. In the case of the W mass, we 
performed fits using KorWan to the invariant-mass distributions from KoralW&YFSWWS 
and RacoonWW, where various effects have been switched on and off. We found, for ex- 
ample, that results of the two programs differ in terms of the fitted by < 3 MeV. 
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We also performed a scale-parameter estimation of the missing effects. From all these 
studies we estimated the theoretical precision of Mw for the above programs at ~ 5 MeV. 



All this is described in detail in Ref. [0. A similar analysis for the TGC A was done 
in Ref. In this case the fits of A to the cosB^y distributions were performed using 
not only KorWan but also dedicated polynomial fitting functions. These fitting functions 
were constructed using MC-genarated data from KoralW&YFSWWS and RacoonWW, not only 
at the parton level but also including the ALEPH detector-simulation effects. Based on 
these studies, we estimated the theoretical precision for A at ~ 0.005. 

6 Conclusions 

There are two independent Monte Carlo programs for the precision predictions of Im- 
pair production at LEP2: RacoonWW |0|-|3 and KoralW&YFSWWS |7|,PH2^. They include 



the 4/-background contributions as well as all the necessary radiative corrections at the 
precision level required by LEP2. The agreement between these programs for the main 
observables is within the required accuracy of the LEP2 experiments. From comparisons 
between these programs, comparisons with other calculations, and investigations of vari- 
ous effects, we have estimated the theoretical precision {6th) for three of the main LEP2 
WW observables: Sth^ww ~ 0.5%, SthM^ ~ 5 MeV and 6th^ ~ 0.005. Comparing this 
with Table 1, it is clear that this satisfies the LEP2 precision requirements. However, 
for the future linear colliders (LC) this is not sufficient and further improvements in the 
theoretical calculations are necessary. 
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